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ABSTRACT 

This thesis is concerned with the steady state response to a step 
input of a non-linear servomechanism having viscous friction and a load 
position unity feedback loop. A gear train with backlash between re- 
silient gear teeth was located between the motor and load. 

The physical equations describing the system are presented and 
adapted to a form for the Control Data Corporation 1604 digital com- 
purer. Theseemputéreprepran vecd im this amalysis is préserited along 
with its flow diagram and description of its operation. 

The results’ ofthis" thesis are presented in several forms using 
the»parameters of system damping coefficient, distribution of friction 
and inertia, backlash angle and coefficient of restitution. The re- 
sults are applied to sample problems in illustration of their use for 
design and analysis. 

The authors wish to express their appreciation to Dr. George J. 
Tiweiter of the Department of Eleetriecal Engineering, and to Dr. William 
Wainwright and Mr. Edward N. Ward of the Department of Mathematics for 


their assistance and encouragement in completing this work. 
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TABLE OF SYMBOLS 


Se (ONE )* 

an (THETAL ) 
©. (THETA DL) 
E 

K (KMOTCONST ) 
Tm 


(RHO ) 


(MOINERT ) 


(FMOTOR ) 


q, 

Tum 

e 

Jim 

om ( JLOADVAL ) 
J 

Fim 

st (FLOAD) 
A 


(DELTA ) 


Angular step input (radians). 
Desired load displacement at steady 
state 
Output displacement of load (radians). 
Output velocity of load (cadens) 

sec 
Error in radians. E- 0, - GQ. 
Constant of proportionality relating 
position error to torque developed by 
the motor. 
Developed torque of motor. 
Torque of load. 


Tongue of lead meferpred tometer. 


Number of teeth on gear 1 
Number of teeth on gear 2 


Inertia of motor measured at motor. 
Inertia of load measured at load. 


Total inertial measured at motor 
J = Jm +e° dL 


Friction of motor measured at motor. 
Friction of load measured at Load. 


Backlash of gearbox measured in radians 
at the output. 


Slope of phase trajectory (combined 


system). a Gee _ i. 
MMs * Fae 


Slope of phase trajectory (load floating 
free). jie d be. = 4 


ee a ee 
oe 


Ia. a 


c 
Velocities of motor and load after impact. 


*( ) Terms indicate computer mnemonics referred to in Sec. 3, Computer 


Program Development. 
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(RESTITUT ) 
( ZETA ) 
(OMEGAN ) 

( OMEGANSQ ) 


(FLFTPRIN ) 


(INERTRAT ) 


(THETAM ) 


(THETA DM ) 


Coefficient of restitution. 
System damping coefficient. 
System natural frequency. 

System natural frequency, squared. 


Load friction ratio with respect 
to friction of the systen. 


Inertia ratio. 


Displacement of motor (radians). 


Velocity of motor —- 
sec 





1. Introduction and Background 

The problem under investigation was the response to a unit step 
input of a nonlinear servomechanism having viscous friction and load- 
position unity feedback loop. A gear train with backlash between the 
resilient gear teeth was located between the motor and load. A de- 
termination of the possible existence of a limit cycle with variation 
of several parameters was made. The variable parameters considered 
were motor, load, and system time constants, backlash and the co- 
efficient of restitution of the gear teeth. If a limit cycle existed, 

a further investigation was made into the change of the size of the 
limit cycle with variation of the above parameters. 

If a gear train is required between the motor and load ina 
servomechanism, the gear train may be treated as ideal with perfect 
meshing of the gears. In this case the components of the system are 
joined at all times and the system may be described by a single linear 
Girercnitial equation. “The classical Linear solution in”eruber the 
frequency or time domain is the result. 

Present production methods cannot meet the requirements of an ideal 
gear train, having perfect engagement of the teeth, nor is the ideal 
gear train desirable from the standpoint of wear on the gear faces. 

The result is that the practical gear train has separation of the gear 
teeth when the velocities of the motor and load are different. The re- 
sponse of a servomechanism with backlash is in effect discontinuous or 
nonlinear. The net system acts as three individual but dependent systems 
with the boundary conditions of each being specified by certain physical 
laws. 


Two accepted methods of analysis are available for solution of the 
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problem. The first is the solution by describing function methods 
while the second is by analytical methods. Describing function 
methods are primarily frequency response techniques and depend upon 
the assumption that the nonlinear element may be considered as linear 
over the range of consideration. A further assumption is made that 
the input to the nonlinear element is a pure sinusoid with no harmonic 
frequencies. The describing function method is primarily concerned 
with steady state results. 

The-analytical solution, which ds thé second method of considering 
the problem is primarily a transient response technique concerned with 
solution in the time domain. This approach also requires certain 
assumptions concerning the behavior of the system. The assumptions 
used in several investigations are pointed out in the following para- 
graphs, while the assumptions of this investigation are outlined in 
section 3, Development of Equations. The validity of the analysis by 
any method depends on the accuracy with which the real system has been 
defined. 

Chestnut and Mayer, Ref. a, have considered the backlash problem 
using describing functions and supported the findings by analog com- 
puter studies. The main difference between this cited work and previous 
describing function analyses was the simulation of backlash by a dead 
zone between the motor and load and springiness in the interconnecting 
shaft. 

Various analytical approaches have been made to the problem of 
backlash, Lutkenhouse, Ref. b, used graphical methods in the phase 
plane to solve several cases of plastic impact between the gear teeth 


of a second order servomechanism. Pastel and Thaler, Ref. c, developed 


analytical equations for the existence of limit cycles using phase 
plane equations for the caseof plastic impact and no load inertia. 
Knoll and Narud, Ref. d, investigated limit cycles in the phase plane 
using an analog computer. The investigation of Knoll and Narud 
covered a wide range of parameters for the case of plastic contact 
between the gear teeth. 

To analytically describe the total system response of a second 
order servo with backlash, three differential equations and one or 
more algebraic equations are required. The system as a whole may be 
treated as being piece wisé linear. One differential equation is used 
to describe the entire system when the gear teeth are in contact and 
the motor is driving or braking. When the gear teeth are not in con- 
tact, two more differential equations are required, one for the motor 
alone and one for the load alone. The boundary conditions of the dif- 
ferential equations are determined from the solution of one or more 
algebraic equations, expressing the laws of conservation of momentum 
and energy. For plastic impact with no bounce of the gear teeth, the 
law of conservation of momentum must be satisfied. For perfect elastic 
contact between the gear teeth the law of conservation of momentum and 
the law of conservation of energy must be satisfied simultaniously. 
Intermediate cases between perfect elastic and perfect plastic contact 
can fulfill only the law of conservation of momentum. However, total 
energy accounting may be made in the intermediate and perfect plastic 
impact cases. 

New, Ref. e, adapted the differential equations of system motion 
and the law of conservation of momentum (plastic contact) to a digital 


computer analysis of a second order servomechanism with backlash. The 





digital computer was chosen as the method of solution for this investi- 
gation primarily because of the flexibility of data presentation and the 
ease with which cyclic parameter variation could be obtained. 

In later sections of this thesis, the differential equations for 
the system and the algebraic equations for the impact boundary conditions 
are developed. The computer program and its flow diagram and modes of 
operation are pointed out. The results of the solutions of problems 
are presented in Several fomms and observations are mede in Seetien 5, 
Results and Discussion. Sample applications of the results to mechani- 
cal systems are presented fn Section 6, Application of Results. 

At this writing, associated work in the area of transient response, 
primarily peak o-ershoot and settling time, of a second order servo- 
mechanism with backlash and plastic and elastic impact is being pre- 
pared as a thesis by C. E. Andrews and R. A. Kelley at the U. S. Naval 


Postgraduate School. 





2. Development of Equations 


The equations developed were those of a phase plane analysis. 


Prior to defining the net system equations, the assumptions under which 


the analysis was made will be stated. When necessary these assumptions 


Will be ampiified and referred to later in this work. 


It was assumed that: 


ue 


The gear teeth were initially in contact and the initial 
conditions of the system were all equal to zero. This was 
later proved to be an unnecessary limitation for the study 

of steady state reSponse. 

Plastic deformation of the gear teeth during steady contact 
and impact and any torsional deformations of driving shafts 
are negligible. 

The inertias of the gears and drive shafts are considered 

as part of the load or motor inertia depending on their 
attachment in the system. 

The law of conservation of energy was completely satisfied 

in only the perfect elastic case by maintaining the total 
mechanical rotational energy of the system constant at the 
instant prior to and after impact. The law of conservation 
of mementum is satisfied in plastic, elastic, and intermediate 
cases. When the law of conservation of momentum is the only 
equation required to be satisfied, the energy lost from the 
system is dissipated in the heat of infinitesimal deformations 
of the gear teeth. 

The gear teeth are in contact only instantaneously during 


impact for the elastic and intermediate cases (excluding 








plastic) and that the impulse torques of drive, friction and 


bearing supports, etc., are zero during impact, e.g. 


\ a. dt =O and ©O= 6’ 


6. The coefficient of restitution of the two opposing gear 
teeth is the same or is described by an equivalent co- 
efficient if the two gear teeth are of unequal coefficients. 

7. Backlash is assumed to be equal at all points on the gear 
circumference. Backlash is measured at the output shaft. 


A block diagram of the system considered is presented below. 
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The equations for the wi and load in contact are: 
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The equation for load alone with gears nee in contact is: 


Ju Oc + hu be =Q 
or a s. .. as 
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The equation for the motor alone with the gears not in contact is. 
Rae Sy ie FL. 8 ara = K (6 Be. 6.) 


The method of analysis used by New, Ref. e, is used in order to 


be able to examine the system independently of its Con 
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and making the indicated substitution to a transformed (*) coordinate 


system for the system equation 
* 2 ger 2 A 
= Wn Oe 


wer Ge + atin Ge. + Wn Oe 
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and introducting the slope equations of the phase — 
© R results. 
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If the load equation is first put in thé form for the phase plane 
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~ 
O results. 
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The same result may be obtained by setting ()), =1. The results of such 
a transformation require inverse scaling for practical application, 
examples of which given in Section 6. It is noted that the slope of 
the load-free equation is mathematically the same whether in the phase 
plane or in the transformed phase plane. It is pointed out at this 
time that the equations later developed to satisfay the laws of con- 
servation of momentum and energy are independent of the system natural 
frequency. 

To establish equations for the law of the conservation of momentum 
and a relationship satisfying the law of the conservation of energy, 
the following schematic is used. 
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By representing the inertia of the load and motor as inertia of the 


gears the following figure is obtained at the instant of impact. 
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The torque equation may be expressed ibe by 6 a a de. 


Impulse is equal to the time rate of change of momentum. The ex- 
pressions for the rates of change of momentum of the gears treated 


separately may be written: 


of 
) | “Om é 
we 7 0 ec CO 


—T, (iF ildt = Vm ( O., - b.,) 
At >o 
= ),\Frlat = Sm 8m - Sm) 
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where the primed velocities are those following impact. 
If the momentum is to be conserved in the system the impulse 


functions are equal and opposite. 


in [8m - Om) = ota (8) - 6) 


ma Tm 1 ( 6 - Om) = (be - 6) 
Siegal? 


from substitution of the equations i 
6 
} ae 
Cr Se = e 


The assumption that all other impulse functions are zero at the instant 
of impact is restated at this point. 
A definition of the coefficient of restitution will be made for 
purposes of this work: 
. » | ci] an 
ee coe ee © bimn- be) 
(Su + 8, \ CP om - 6.) 
It will be shown that this definition of e will satisfy the law 
of conservation of energy. Consider the case of e=0, , 
e/ / a.) 
b=-(POm- Se) oo = oO, 
in which case for plastic impact (e=0) the gears are moving at the 
same velocity following impact. 
The expression for the coefficient of restitution is examined for 
the case of perfect elastic contact, e=1, for which the law of conser- 
vation of energy is satisfied. 


Solution of the restitution equation for e=l yields 


1 








PC Sn- bc =- p bms Oe 
e (Om + 8 wn) = yi + Be 
For conservation of —s : a P 
L (3m mt Tbe) = L(TmOqrd 6c) 
2 TS | : 
. Jin (8m = Om) = (8? - 22) 


and factoring 
(iit Sy) (4° bm) = ( det 6) [bc 4) 
The momentum ae which is independent of ¢ is: A! 
be: (et 6.- © e) 
or rearranging e + len ( B yy7 6) ( 4/- 6 .) 


The factored energy equation ‘r oe divided by - momentum 
equation to yield (( Sivagat 8) = (6, a .) 
which checks with the definition of the coefficient of restitution 
for e=l. 

The momentum equation implied from the impulse approach that 

d+s> 0 hence 6.= 6. and 6 = 6) . The definition of the 
coefficient of restitution is independent of position. It may be 
noted that the equations for energy and momentum depend only on the 
inertia distribution, gear ratio and angular velocities. 

In manipulation of these equations it is seen that load inertia 
may be transferred to the 6 » shaft by multiplication e* » the same 
transfer may be accomplished with load friction. Since 6, » ome, 6, - 


are related to oT and 6, it may be reasoned that the results 
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ror C = 1 may be extrapolated to physical systems where (= # ile. 
Since a total energy accounting may be made for the system with 
the equations used, the rotational energy lost is attributed to the 
heat of deformation of the gear teeth. This transfer of energy could 
in fact be determined from the equations used. Thus the conservation 
of rotational energy e=l is a special case of a broad interpretation 


of the law of conservation of energy. 
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4. Computer Program Development 

The physical equations of the net systems were programmed for 
solution using the Control Data Corporation 1604 high speed digital 
computer utilizing paper tape program input and magnetic tape output. 
An IBM 717 line printer was used to extract data from the magnetic 
tape. The Control Data Corporation machine library and the U. S. 
Naval Postgraduate School computer subroutine library were used for 
assembly, Runge-Kutta-Gill numerical integration, and decimal output. 

Several changes were made in the forms of the physical equations 
of the net system in order ae duplicate computing operations 
and to fit the equations to a form suited to the variable parameters. 
In the table of symbols, computer mnemonic (m ) terms which are used 
in the assembly subroutine have been indicated by parenthises and will 
be defined when encountered. 

The equation of the system with motor and load combined 

OW. = wey, Qo ’ ay Ce. —~ A S C1) v1 Oe 

was put in the form 


oe Ldn Be = Oy Oe. ~# 6. Where 
A- 25 Wn and 


THETAL = Oe WE = Oe 
THETALD = \\ = Oo ae 6S 
vor = Oy ommcansa = (4) 


were the computer m terms used. Underlining will be used in this 
section to denote the transition to the m terms of the camputer 


equations. 


The differential equations were solved using the Runge-Kutta-Gill, 


7 





(RUNGE), numerical integration method. The RUNGE subroutine re- 
quired the definition of synonymous terms for the four iterative 
cycles used to produce one extrapolated set of variables. The incre- 
ment of the independent time variable chosen was 0.01 sec. for all 


> O.| amd 0.007 sec. for all Ded Oe 


The load free equation 


Se i sane 1% a ls, was put in the form 
oe Be : 
or 4. B S C. = O for the computer. 


A separate set of RUNGE synonyms was used for the load free equation 
in order not to destroy the ee, computations before they were 
determined of no further value in computing and also in the making 
of program decisions. 
O, was designated V and On was designated VDOT for the computer 


program. 


The equation for the motor alone 


+ G. 6,,+k & =k oe 


oa 


was put in the form 


a4 - 


and for RUNGE, the terms 


©», = THETAM Oy), = THETADM = W 
- 8 
and Ow = WOT were used. 





The equation representing the law of conservation of momentum and 


the definition of coefficient of restitution were combined to the forms 
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Bo tate: e 6,,, (1+ €)+ o(Pb-9 
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rete CPUmMt eC Ge 
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for the computer program. No provision for additional m terms was 


o| 
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made to denote the primed values. The additional terms: 


e=RESTITUT P= RHO Gnd P= RHOSQ 


~ 
* 


were used. 

Iwo equations were used to define the boundaries of operation of 
the load and motor when they were acting separately, 

Sg eC Sm Ond B.= Pp SmtA 

Two major decisions of the computing cycle were the determination of 
the point where the motor and load would float free, FLOATEST, and 
the response of the system to impact of the load and motor gears when 
the boundary conditions on position were met, COMBTEST. 

The first decision, the point of float free, is made on the 
equality or inequality of slopes in the phase plane, i.e., at separation 


N= No. Equating the slopes: 
ie 
Wy (bg- Oe) -4 ri Cy = © 
me. 
The above equation is solved after computing each point in the 


combined phase trajectory. If NV a > O the system re- 
mains combined. If the system separates, havin 

Nis +N, < © the system separates, having 
the response motor alone and load alone immediately thereafter. 


The second major decision, that of whether after impact the 


LY 








response should be that of the combined system or the response of the 
load and motor acting separately, is made on the basis of the resultant 
velocities after impact and upon which side of the backlash the motor 
and load positions are found. 

The computer program, its flow diagram, a description of the pro- 
gram operation, and an explanation of the readout of information with 


examples are presented in the Appendix. 
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4, Methods and Scope of Examination 

Iwo modes of operation of the computer program were used. The 
first mode considered was the phase trajectory. This mode was used 
to check the operation of the computing cycle. Time, velocity and 
position of the load were printed when the system was combined. When 
the system was separated, the values of motor velocity and motor 
position were also printed. Although the time increment used for 
computations was either 0.01 sec. or 0.004 sec., printed outputs 
for the first mode were taken only every 0.1 sec. Additional printed 
outputs were taken at the time of contact of the gears just prior to 
and immediately following the solution of the momentum and restitution 
equations. A typical printout of the phase trajectory mode of operation 
is shown in Appendix C. 

The second mode of operation was used for printing only the maxi- 
mum computed positive load overshoot position for each cycle, the 
associated load velocity and the exact problem time of the computation. 
A sufficient number of print-outs for each problem was obtained to 
prove limit cycle existence and average size. Since this investigation 
was concerned primarily with the steady state response, the second mode 
of operation was the one utilized to obtain the majority of the data. 
This read-out method markedly decreased the data reduction time for the 
problem analysis, as opposed to the analysis of steady state conditions 
provided by a full phase trajectory print-out. Since the time required 
for computer read-out was several magnitudes greater than the computing 
time, valuable computer time was saved by this mode of operation. An 
example of the printed output obtained when only maximum load position 


overshoot was of interest is given in Appendix BP. Approximately 1400 
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phase trajectories were solved using this mode of operation. Each 
solution required an average of three minutes of computer time. 

For the purpose of general examination of the problem it was 
assumed eer? =1,K =1, Os 1, Sat Desie = 1. Specific solutions 
were made with C 2 ol , Uj# 1 to determine scaling effects. The re- 
sults of these solutions are presented in Section 6, Application of 
Results. 

The program can also be made to solve linear systems by setting 

AERO: 
The major variable parameters used in this investigation were: 
5 NOL 14 OL2 003, OCS. 5556. ee a 


i“, O , O82, Hh, OM 0.8, 10 
a 


ran 

_St= 0.1, 0.2, 0.5, 0.8, 0.9 
Im O.9, ES Dd, Oeey. ec. 
A = O.4 radians 


eure 0-0; 070,.1.0 
Additional parameters of 

Cf - 0.04, 0.1, 0.9; and A = 0.01, 0.03, 0.01, 0.15 radians 
were beat in certain instances to examine particular characteristics of 
the response. 

The program was limited somewhat in that the value of an = 50 
and zero were exluded due to generation of undefined a Es 
quantities. 

The value of backlash was made abnormally large, A = 0.3 radians, 


to allow an easier interpretation of the non-linear response. Since 


the influence of A was linear, this caused no inaccuracies. All 
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graphs for limit cycle size were plotted with /\ = 0.3 radians. 
The effect of values for /\ other than 0.3 is also shown in graphi- 


cal form. 


ae 








5. Results and Discussion 

Two types of system response were observed. The first type was 
a convergence to zero error, which was a case of static and dynamic 
equilibrium of the system. This is a characteristic of a linear 
second order servomechanism. The second type of system response 
observed was the divergence or convergence to a stateof dynamic equi- 
librium. This second type of response is termed a limit cycle and is 
characterized by cyclic travel of the load through the same points 
of the phase trajectory. When the second mode of camputer operation 
was used, both the maximum load position per cycle and the period 
of the recorded position were analyzed to determine if the character- 
istic response was a limit cycle. Exact repetition of cyclic values 
was impossible due to the numerical methods used and the computer 
round-off error for the output routine. Typical examples of the 
limit cycle and no limit cycle response are given in Appendices 
D1and D2, 

No case of dynamic instability or divergence of the system with- 
out bound was observed in the investigation. Intuition might lead 
to this same conclusion when the following unique characteristics of 
this system are pointed out: (a) Asa maximum limit, energy was 
conserved at the instant of impact. (b) Energy was supplied to the 
system only in finite amounts and at a finite rate by the error 
detector. (c) Energy sinks were present in the system frictional 
elements while energy storage units were present in both the motor 
and load. (d) No undefined limits appeared in the applied physical 
equations, 


(es 
Au = © and oD excluded. 


Ce 


The results of the investigation are offered in the form of 
charts for parameter areas of limit cycle existence and nonexistence 
im Fig. 1 for ¢ =s0., 0.6, Ov8 eand)Fie. 2 forves= 1.0. The para- 


2’ 
meters ( im and p Ie are designated abscissa and ordinate 
ray ym 


respectively. The parameter points of examination for which solutions 
were Obtained are circled. The zone enclosed by a particular a line 
is largely an interpretation by the authors of the results of the in- 
vestigation, and shows the area in which no limit cycle existed. Due 
to the large change of the variable parameter values between solutions, 
the charts are not exact. All circled intersections interior to the 


&, ) 
. line indicate parameters which resulted in no limit cycles, e.g., 


im. 1 for e = 0, 0.6, 0.8 at & pS, = 0.4, the values of Pee = 
TT 
C6, 00S at a p° Te = 9,0 and the value of Naa =O ./even & 
ym or 
oO Nu = 4.0 resulted in no limit cycle. Similarly, for “§ = 0.5, 
Jm 
the values of Pfu = 0.4, 0.6, anda ere = 1.0 as well as all 
er ym 
other parameter intersections interior to the S = 0.5 curves, re- 


sulted in no limit cycles. 
All circled parameter intersections exterior to particular Ss lines 
define paramters which were examined and resulted in a steady state 
limit cycle. ‘The area to the right of the single \ = 1.0 lime 
designates the parameters for which no limit cycles were observed. 
Limit cycles were observed for all ine 0.3 for the parameters examined. 
The investigation was made with A = 0.3 radians; however, re- 


sults of the investigation of variation of limit cycle size indicate 


ao 





that existence of the limit cycle is independent of the size of back- 
lash in a nonlinear system for a A greater than zero. Variation of 
A had the effect of a linear variation of limit cycle size only. 

The zones depicted on the existence charts are approximately 
represented by straight lines; however there is reason to believe that 
a complete examination would show a curvature in the lines. This 
observation will be elaborated upon later when the relationship between 
the existence charts and the figures of limit cycle size is discussed. 

From Fig. 1 it is seen that there is only a very small difference 
in the limit cycle existence zones for plastic impact (e = 0) and the 
intermediate values of the coefficient of restitution, e = 0.6 and 0.8. 
When comparing the elastic impact (e = 1) in Fig. 2 to the plastic and 
intermediate cases in Fig. 1, the most obvious point of camparison is 
the drastic reduction of the areas in which limit cycles did not result 
for e = 1.0. However, it may be generally observed that with any e 
held constant, the area in which a limit cycle could not be obtained 
increased with the increase of systen i : 

For all values of e, there is an apparent common enclosed zone 
centroid which is described by the system parameters Cr F, mO.5 

a 
and . J, _- 


stants iia = Fn/ dm5 since OU. = 1 and ( = 1 were chosen for the 


1.0. These parameters result in the system time con- 


investigation. There is also an apparent symmetry about the centroid 
of opposite quadrants which can be related to the fact that since 

UJ, 1 and Ps 1, the points diametrically opposite across the 
centroid of the figure (using the decimal values of the parameter inter- 
section) represent an exact exchange of component time constants. (See 


Fig. 2 for example). Two diametrically opposite points may have 
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exactly the same motor and load time constants; however, they are 
descriptive of two different systems since the friction and inertia 
distributions are different for the two points. 

Figs. 3 through 26 present the maximum positive error in radians 
of the limit cycle for a unit step input and with backlash equal to 
0.3 radians. The use of the charts is explained in — 6, Appli- 
cation of Results. Each figure is presented with er L/*y and e 
constant. The independent parameter er) J is indicated as abscissa 
and curves of sonasm ‘5 are plotted with limit cycle error as the 
ordinate. The results were presented on a log log plot since the 
observed limit cycle errors encompassed a feage from “0.01 to 170 
radians, and the appearance and disappearance of limit cycles and 
the change of limit cycle error was more readily apparent on the log 
Wee form. 

In @eweral as e increases, the size of limit cycle decreases, all 
other parameters remaining constant. However this is not without a 
few exceptions * the generalization is most accurate for large size limit 
eycles. It may be noted that when viewing the individual charts, if a 
limit cycle exists at very low C5: /s, ratio there is a tendency for 
the same system to have a smaller limit cycle for higher @ a,/3, Pave. 
In addition, for these same conditions, if limit cycles exist for 
several values of S » there is a convergence toward the same size of 
limit cycles for each S as Og, /5, increases. At the value 

Ot, /Fy = QO, a limit cycle will always occur for any value of e or 
S. . 
Viewing the charts in sequence, with either e constant or ‘entian 


constant, it is noted that there is continuity of pattern flow between 
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adjacent charts. This continuity between charts appears to form a 
symmetrical pattern of constant ‘S lines which can be interpolated for 
intermediate values. When the figures for Oe /r, =—wOwie emideO26 
are analyzed, this symmetrical pattern is observed in each chart about 
Os /5 = 1. The limit cycle decreases and disappears at the 
opposite extremes of p3,/z.. The charts of pt,/F,, = 0.4 and 0.6 
are symmetrical to each other, while e £,/F rp ="Ov2 andeOvws ame 
similar. At the extreme values of yi £ Py = 0 and 1.0, no symmetrical 
comparison can be made. An additional point of symmetry about 
e t,/F, = 0.5 can be visualized for all values of e. These para- 
meters for symmetry are noted in the discussion of the existence charts. 
The irregularities in the curves of medium values of S can not be 
explained by the authors; however, it may be mentioned again that the 
points of symmetry are characterized by exact exchange of component 
time constants. 

The noted symmetry of the limit cycle size charts shows a direct 
relationship to the symmetry found in the existence charts. Since the 
existence charts show contours of constant S for a stirface,of zero 
limit cycle size, one is led to the possibility of a three dimensional 
presentation to show a variation of limit cycle size for the same para- 
meters used in the existence charts. Such a three dimensional figure 
for constant e using e f/F PI/5, and. 5 as the axes is shown 
in Fig. 27. Other surfaces could be added above this zero limit cycle 
"Ploor" to show the variation of a given limit cycle size by contours 
of iS . The ultimate result of this line of investigation would be 
the evaluation of the poles and zeros for a conformal mapping presen- 


tation. Unfortunately, time did not permit the authors to investigate 


26 


this avenue further. 

The results of the investigation of the dependence of limit cycle 
size on the amount of backlash are presented in Figs. 28 and 29 using 
coefficients of restitution of O and 1.0 ata S = 0.6. A thorough 
investigation was made at this S Since it is an average system 
parameter often used in design work. The limit cycle size was found 
to vary directly with backlash size without exception between the 
values of 0.4 and 0.01 radians. When limit cycles existed at A ="Oos, 
there was no disappearance of limit cycle with decreasing magnitude of 
AA. Where no limit cycle existed, at /\ = 0.3, limit cycles con- 
tinued to be non-existent throughout the range of A examined. 

Several solutions of problems with < = 0.6 and e = 0.6 &0.8 and Ss = 
O.1 &e =1.0 & O with various /\ values, proved the linear variation 
of limit cycle size with backlash size for all values of e. 

When a servomechanism under steady state conditions is disturbed 
by a small perturbation, it will return to the original steady state 
conditions. However, the transient response due to the perturbation 
may be quite different from the transient response which occured due to 
the original signal input. The concepts of energy and static and dy- 
namic equilibrium mentioned in Section 1 may be better understood as a 
result of the perturbation technique if the usual initial conditions of 

8. (0) = § (0) = 0 are considered as a disturbance fram the stability 
conditions predetermined by selection of system parameters and signal 
imput. 

This same small perturbation technique was simulated mathematically 
using the digital camputer to obtain several solutions. The cases 


examined and their initial conditions were: 


oa | 


o-d 
p Q 
™, icecrrrmr 
© Oo 
i 
1] | 
© © 
tl 
¢ © 
H oo 8 
re i=-—™ 
© Oo 
Ne Nie 
{| il 
a 
© NO 


system separated. 


ae rs 
b) = 0.6, PI /J, =0.1/0.9, Pr /F,=1.0,e=1.0, A = 0.3 
Oe (0)=1.0 , Ge (0) =1.0 
Om (0) =1.0 , @ m(0) =1.0 
system combined 
> ei 
c) = 0.4, eee = 0.5/0.5, Of ,/F =0.4, e=0.6, A=0.3 


S 
Oe (0) =o, Oc (0) =1.005 
G 


m (0) 0 Gm = 1.00 


system separated 

In the three cases examined, the same limit cycle was obtained, 

both in size and period, as that resulting from the normally used 

e 
initial conditions, 6 e (0) = 6 ec (0) = 0. The maximm difference 
between the perturbation response and the original limit cycle size 
for the three cases was 18 x 10°” radians. Two types of transient 
response were observed. These were convergence to limit cycle in the 
first two cases and divergence to limit cycle in the last case. 

The results obtained by the small perturbation technique supported 
the validity of the analysis of the problem by digital methods. This 
result also concurs with the findings for the case of e = 0 by Knoll 
and Narud, Ref. d, in that the steady state limit cycle size is campletely 
specified by the choice of system parameters and independent of system 
initial conditions. 

Mention is made here that the computer program could not be run 


normally for the perturbation method. The computer was stopped after 
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the print out of system constants, then the initial conditions of the 


servo system were inserted manually. The program was then restarted 
at the desired point for computing, either as a combined system or 
with the motor and load separated. The program was then allowed to 
operate to the normal completion of the solution. 

The work of Ref. d reported that the system transient response 
for the case e = O arrived at a steady state limit cycle by two 
possible phase trajectories. One of these phase trajectories was 
smooth convergence to the limit cycle, while the second was a converging 
overshoot toward the inside of the limit cycle followed by a divergence 
to the limit cycle. From the results of this investigation for the 
eases e = 0, 0.6, 0.8 and 1.0, similar transient indications were ob- 
Taine . 

Several computer solutions were accomplished for the cases A = 0.3 
mediens, €= 0, 1.0, R + 2.0, Ly) n= )..60,;, Ke= J, # 1.0,. im addition 
various combinations of P # 1.0 and Ly) n # 1.0 were obtained. The 
purpose of these solutions were: 

1. To determine the general applicability of the computer program. 
2. To verify the scaiing developed in Section 2, Development of Equations. 
4. To determine the applicability of the results of this thesis. 
The following results were determined: 

A. Two solutions were obtained (e = 0 & 1.0) for the case < =m sf ey 

@= 1-0 Lb) ® =4..0, mS 2.05 = 2G tu = 0.8, P- 4.0. 

a om 
The results of these solutions were the same (time, load velocity, load 
position, and limit cycle size) as the solutions obtained for at = One 
a 


ed 
Pimalie = 2 10, < = 0.4 using the nondimensional parameters br) he= 170; 


Jb ™M 
= er JI = 1530. K lee). 


2? 








B. Solutions were obtained for e = 0& 1.0 for the case as mEOe, 
2 
= ar, Qn —ETO es = OK Seay P fe =a.) ye = m0. 
FT ym 
The results of these solutions yielded the same size maximum overshoots 
and limit cycle size as the nondimensional solutions for Ss = 174, 
Ve. oe 
a saOrlp, R yt=4.0. The times of maximum overshoots in 
wi Tm 
this case were one-half the times of occurence of the same overshoots 
for the nondimensional solution. The velocity at a maximum overshoot 
was twice the value of velocity obtained for the same overshoot in the 
nondimensional solution. This same comparison can be made for Cases A 
and B. The results of these and other solutions prove the validity 
formulae used in Section 2, Development of Equations, concerning the 
transformation to the coordinate system. (*) 
C. Solutions were obtained for the case e =O0& 1.0, as me Ones, i? = Oe, 
2 
ba) n= EO, a7 = On8:, Ju = LyeO% care = Oe2, © Je — Lae. The 
Ee Pr tm 
results of these solutions were the same (time, load velocity, maximm 
overshoots and limit cycle size) as for the nondimensional solutions 
2 ra 
a =0.4, e =1.0, G)n=1.0, Pte = 0.2, PJ. = 1.0. 
eT Tm 


D. Solutions were obtained for the case e = 0 & 1.0, Ss = (Oh 


et, LO) B= 250; Fu = 0.8, J amo, oo =0.2, PrIe= 1.0. 
aa Im lea 
The results of these solutions yielded the same size iron overshoots 
and limit cycle size as the nondimensional solutions for ss ="Onas 
Pl. One, e° x = 1.0. The times and velocities in this case were 
i 
scaled by W n as rs Case B. 
Several phase trajectories using e = 0 were compared with those 
obtained by New, Ref. e. Generally good correlation was obtained on 


the value of limit cycle error and in the transient response. Mall 


disagreements were expected, resulting primarily from differences in the 
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programming procedures and analysis of results by read-out methods. 

Sample phase trajectories for the cases of e = 0 & 1.0 are pre- 
sented in Figs. 530 and 31. Although the transient analysis of the 
system is beyond the scope of this work, these figures may aid in 
understanding the physical problem. 

Of academic interest was the accidental operation of the computer 
program for several solutions using S = 0. Several responses were 
observed. One was a nearly pure oscillatory system from the first 
overshoot and thereafter another was a very slow convergence or slow 
divergence of the oscillations after the first overshoot. The con- 
vergence or divergence probably resulted from errors generated in 


numerical iterations and approximations. 
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6. Application of Results 
In this section, the results of this thesis will be applied to 
hypothetical physical systems. Methods will be indicated by which physical 
systems may be analyzed or designed. Attention is called to the assump- 
tions stated in Section 2 under which this study was made. 
it is desirable to find an equation relating certain system 


variables, therefore the equations: 





Ymt pd am | 


7 ‘wane =: +H 
(2) AGidy = Tete te = seal hed 
Tint Se —_s p Je 


will be manipulated to the form 


and 





=” 
+b: pL iw ‘a 


(3) ~ : 
Seely + eae] Ve VR “ 


a: Gm 4) 


= 
i § — | 
oF. Tm & (1+ ¢ Se ) 


In view of the form of Kquation (4) Tor \ and recalling that the 


if 


2 
variable parameters used on the figures were R a and \ f. = 


a oe wm 
L 


= oe ra it is seen that the expression for a. 
fant PEL TEL 


in fact includes on other two variable parameters. 


B= 





From discussion concerning the effects system natural frequency, 
it will be recalled that the value of overshoot and hence limit cycle 
Paesnot affected py (A) n; however, the time and velocity are scaled 
for Ww n#1.0. Thus the transient response characteristics may 
be related to the same parameters used in these applications. 

If a motor and load are selected for a system and the limitation 
is imposed that system operation shall not result in a limit cycle, 
Equation (4) and the limit cycle existence charts may be used to 
determine the unspecified system parameters x » ee and % . The 
coefficient of restitution for the proposed gear train material may 
be determined by the use of the equation in Section 2. 

To proceed, select a value of K and e , then determine the value 
of 3 from Equation (4), Enter the limit cycle existence chart for 
eae value of e, system < and the parameters po Se and a ° 

Sm alae 420) 

ZF 
Determine if the selected parameters describe a point of no limit 
cycle. Since limit cycle existence is independent of backlash (A ra 0 
excluded), backlash is not a variable parameter. Successive trial 
values of the unspecified parameters may be required to meet the 
limitation of no limit cycle. If it is also required that certain 
transient response characteristics be met, these conditions could be 
examined by a similar method. 

In the event that no acceptable parameters are found to satisfy 
the system requirements for no limit cycle, the figures for limit cycle 
size may be examined to obtain a minimum size limit cycle. The method 
of solution is similar to that previously deseribed. 


Example: The material intended for use in the gears has an e = 0.4. 


oe 





The motor and load have the parameters fm = 0.64, Jm = 0.25, f. = 0.64, 


L 


Jy = 1.0 and it desired that the system operate with no limit cycle. 


2 
If the values are selected K = 4.0, eC = 0.5, the values Q Je = 1.0 


and sia ae lieiO' codleeereet Que: O.2 result. saa 
PF eL Fin +. | 
ex are oe. 
Solve Equation (4) for S = .5(.64)(5+1) : , thus 
PastdUit) 7% 


aS = 0.4, Enter the chart for limit cycle existenes, Figs ly 


2 
e= 0.6, \ = 0.4, ose bea = 0.2 C Ji =1.0 


) 
3 ia Im 


It can be seen that limit cycle will exist for these parameters. 








Since K is most easily varied in equation (4), a new K is selected 


at 1.0 in an effort to find parameters for no limit cycle. It is seen 





that 
e z § yt = 04 (2 
eee K= 4,0 
with ail other parameters constant, thus = = hate) Enter Fig. 1 
with = 0.6 = \ 4 7 
- 0. ae = 1,0) | - 0,2 
me Tm 4) 


erFL 
These parameters describe system operation with no limit cycle and are 
well within the area limits. The transient response of this system 
may be undesirable and various cambinations of ¢ and K would have to 
be tested to satisfy the additional requirements of the problem. It 
will be restated that the existence areas of Figs. 1 and 2 are not exact. 
The figures of limit cycle size may be used in the same manner as 


the existence charts. The same parameters of the first case are 





assumed: K = 4.0, P = 055% pe oe = 1503 ! = 0.2 x = Or 
ey Tn + | 
Cis ae 


val 


) 
Ay 








e = 0.6 ) A = 0.1 radians. Enter Fig. lO with the parameters given. 
The ordinate, Limit Cycle Error ( 0-3 | = 0.19 radiens; isecbraidmed. 
The maximum positive error in the limit cycle = 0.13 (0.1) = 0.043% 


(0.3) 
radians. Enter the same figure for e = 0.6, < = 0.8 , PL = 1.0, 
| Ym 

fv, > 0.2 ) A = 0.5 radians; this is the second case of exami- 
Bre. 

nation for existence. No S = 0.8 curve exists; hence no limit cycle 





exists for this choice of parameters for any nominal value of backlash. 
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~. Conelusions 

Within the limitations of the assumptions made and the scope of 
this thesis, it is concluded that: 
A. Two types of system steady state response result. The type of re- 
sponse is independent of backlash size, system natural frequency and 
initial conditions. The types of response are convergence to signaled 
impubeposition.and convergence to a limit cycle. 
B. If a limit cycle exists in the steady state response of a systen, 
the size of the limit cycle varies directly with the size of backlash, 
all other parameters constant. 
C. The greatest change in system steady state response occurred for 
the case when system rotational energy is conserved (e = 1.0). 
D. A point of symmetry for the type of system steady state response 
exists with equal motor and load time constants. For medium to high 
values of system damping coefficient, symmetry with respect to limit 
cycle size exists for an exchange of motor and load time constants. 
E. The results of this thesis are presented in a form which may aid 


in system design or analysis. 
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8. Recommendations 

From the results of this thesis it can be seen that there remain 
many fertiie areas for further investigation. For studies of this 
type it is strongly recommended that the digital computer be used 
because of its speed and versatility. The limits of investigation 
depend to a great extent on the ingenuity of the programmer. Some 
recommended areas of future investigation are: 

a) More precise definition of the limit cycle existence zones 
and a more thorough coverage of the limit cycle size of similar systems 
having the same type of paranieters. A three-dimensional presentation 
as mentioned in the discussion could be obtained. 

b) Include torsion and other deformations of the mechanical system 
in the physical equations of the system to determine their affect on 
the system response. 

c) An analysis of the limit cycle existence zone and transient 
response of a system to determine if there is a relationship between 
the two. An analytical expression for the possibility of limit cycle 
existence might be developed in terms of system parameters. 

d) Examine the transient and frequency response of the system due 
to various Oe input functions. This would require minor modifications 
to the computation section of the existing computer program. 

e) Since the physically realizable quantities of position, velocity 
and acceleration of both the motor and load are available as computed 
quantities, the study of the compensation problem and its optimization 
are obvious areas for further consideration. 

f) The general study of other higher order nonlinear systems is 


suitable for digital computer programming since most numerical inte- 
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APPENDIX A 
COMPUTATION FLOW D/AGRAM 


CDC /604 D/G/TAL COMPUTER 


START LAZYZETA MANSTART 


to STOPNOW, 88 





COMBINED 


from 67, 68,77 





STOP 
NO YES 
a7 
J . ag 
wo to LIMITCYE 
avixcome 
ir . : es 
| 
2z 
ore ON ] 
BIGZETA l 
is i : = 
A | 
| | 
from STOPNOW, 88 [ i 
GOAROUND | u 
iS | | 
i 
16 | 
| 
INOXVALU 
EXITY We———— ad 
18 _— 
onl 
33] : 
NO YES 19 | | 
i te CONT EST.39 
Y 
&O LINBAR No | 
NO VES . | MOTOR 
| 
ai zs Se i Be 
| 
36 | 
i 
Krom 70,7/ 


(¢4 











NS SS ea ee2ee2e2e e2e4 


trons MOTOR, 38 


COMBTESTY 54 


aN Fm 


» COMPLMNT 


NO 
YES NO 
te K 55! 
VN yess N 





Yes 


FLOT PRIN 


med 


ae oe on aoe I 
o<---=--- —<-----,] TF 6 


FLOTSTO 


¥ 

56 
5 
LAS 


Yes 


g 
Se 2k 
— [ #*) 
3 


lo LIMITCYC 


ES 
te 


to COUNT SEC,29 


ee = 


™m 
3 
a 





pLusTesT] 62 


to 
MINVTEST, 
72 











BOUNC from 76 
: 
i 
Y 64 
{ 
j 
i 
Le -—- ™ 

EXITCOUN § 

VAIS PLUSFORSC 
67 


bs = 68 | 


wo 


pREPLANS) G9 


YES No 


D 
‘¢ 


PLUSCORR 


TO 


COM BINED 
RF 


te LOAD,37 


from 74, 75,76 


b> 
L | 
2 


to STOPNOW, 8B 


C5 





from 59,60 


from R7, 43 







MINUTEST] 72 
LIMITCYC 


ON OFF 


to 67 





to 
PRE PLANS, 
69 
to COMBINED, 24 
from 5,6 
Stop Now { 88 
v. _ Exit 
to GOAROUND, 
89 
STOP Jon 
YES NO 


return to program 


99 


STOP AVE 


om ; 





return to progres 


74 








10. 


ae 


2. 


i: 


Li 


APPENDIX A (continued) 
EXPLANATION OF FLOW DIAGRAM BLOCK SYMBOLS 


Clear fault stop. 

Print last overshoot stored in LIMITBUF (symbol 888 ). 
Print out position in phase trajectory of fault stop 
(symbol 666 ). 

Heceor file. 

Actuate fault stop. 

Load zeta from Index 1 (b, ). 

Princ [ast overshoot stored on LIMITBUF (symbol 888 ) . 
Print position in phase trajectory where it was manually 
stopped (symbol 1961). 

Stop if Zeta a's 
BF Wn 


is zeta = 0.17 


A 


is zeta less than 0.1? 

Set up computation time of 0.01 seconds for problem and 

set up for print at every 0.1 seconds when zeta is greater 
wpen, OLE 

Set up computation time of 0.004 seconds for problem and set 
up for print at every 0.1 seconds when zeta is equal to or 
less than O.1. 

Set up limit cycle print-out routine for 20 overshoots and 
average the last 8. 

Set up number of limit cycle print-outs and averaging routine 


depending upon the size of zeta: 


if, 





beg Start Divide Comp. 


Overshoots Averaging bys Time 

20.8 6 \ PaO 0.01 
0.6 if 5 250 0.01 
0.5 8 5 Bae OgeH 
O.4 10 7, one 0.01 
©25 14 10 LO 010i 
0.2 20 12 CuO Of01 
<0.1 “ity 12 50 0. 004 


*The last row printed has this number printed on the far left 
if the problem has a normal solution. This number indicates 
how many overshoots were averaged. The average value has 1.0 
subtracted from it and is printed in the far right of the same 
last row. 

1». Frage reyeet. 

16. Clear buffers and set up for new problem. 

17. Restitution, e, from Index 3 (b;). 
Backlash, ZA , from Index } (b, ). 


Jp, ftom Index 5 (b,.)- 





& from Index 6 (bg). 


lis \(2F On)(Jr) = fis fy =B 
fw = fia 5 eo =a 


18. Have all values been cycled through? 


19. Set STOPNOW after this run to stop absolutely. 
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20. Backlash, f\ , = 0? 


Cl. Print consecutive number of this parameter run and all the 


parameters used: 


aponrA Im fife, e a /? 


rhe i a a ge 
K Yec fm =D Yy=B Wn = 


ee. Set program so always acts as combined system. 
24. Print consecutive number of this parameter run and selected 


parameters for linear system: 
_ + 2 2 2 
2Ppon=A JmtPh KK = Wy fo 


24. Solve the second order equation by Runge-Kutta-Gill numerical 
integration every 0.01 seconds problem time when zeta is greater 
than 0.1, otherwise every 0.004 seconds. Initial conditions 

@ e@ 
a = leai® = = O and = 1.0 or the pre- 
a oO, ? © ©), <, = = 
viously computed point. Next point is computed by four 


iterations of: 


‘° Z ° 

Oc= Wn (Og~ Oc) - 2 PWn Oc 

25>. Set EXIT from print routine. 

26. Load time, eC. and ©, in print buffer. 

27. Go to LIMITCYC to print only the overshoot. Otherwise print 
the phase trajectory points for each 0.1 seconds up to a maxi- 
mum of 12 cycles. 


= 
28. Store O for Ox and Om 


ai 





29. 
30. 
Dale 
ge. 
DD- 


a4. 


aD. 
46. 


Df. 


2. 
oO. 
Wa. 
42, 
43. 
ye, 
i. 
46, 


(Or- poor -2P wn so fet ‘ 


Is problem time 0.1 sec.? 

Is problem time greater than 0.1 sec.? 

Print the number of point computed, time, © ¢ and Oe 
EXIT (automatically set for desired jump-out. ) 

Store time, ©. ; ©. and ©, for motor and load initial 


conditions. 


onige ie 
Is slope Ny N.? 


= positive value? 


Is ©. : 
Store On. - A itr Om 


o& f, , 6 
Solve r= {'|6.) Por CO» and Or 


by Runge-Kutta-Gill. 


= -(Fog On) + Pn)lOn - 2) 


ior Gv, and. es by Runge-Kutta-Gill 
Ou -? Om = @©F 
an ./? Own A = negative? 


pet EXIT from print routine. 


Solve 


pew HXTT from mrint routine. 
Load time, QO. and Oe. 


Load ©y,, and Om - 


Print number of point computed and stored values. 


EXIT (automatically set for desired jump-out). 
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“7. 


48, 
49, 
50. 
5 ii 
Be. 
op. 
a4. 
DD « 
por 
D{- 
BS. 
39. 
60. 
oi. 
62. 


63. 


64. 
65. 


66. 


Solve 


o7 Jm = : “ii & 
6. * ar, pn Om(ire) + Cc(7 — e) 


of 6’ - e( 9 Smt Gc) 


cof fm ed Kk... - JY 
mn - Be i saith ‘i oo 2 


Set EXIT from print routine. 
eo 
e * » 
Is ©, positive? 


Set up to count one phase trajectory when oO, is negative. 

ioe ee => 1.07 

Have 12 cycles of phase trajectory been completed? 

Arrange for no counting until e. is positive again. 
* 0 

Is “70m aa Oe negative? 

Set up — (7? Om i O<) : 

Ts /2 Gm os Oc greater than preselected € ? 

Set up combined systen. 

Is ©, positive? 

Is S. negative? 

Is @.- | positive? 

Is Ge > } negative? 

Is Qc7- ) negative? 

Cz “9 Om= 0% 

Oe- 2 Ome 07 

Is CO. We Em positive? 

Is Oc i Ormpositive? 


ie 
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68. 


69. 


fo. 


fel 


We. 
fit 
Th. 
75. 
76. 
{f- 


78. 
19. 
80. 
on. 
82. 
S5. 


Ont 
86. 
S77. 


Is Slope Ny = No “ 


2 
(1 O-) Un ~ 2 ii, {, 
2 i i ° : 
© = positive value? 

= Lu 
Modify FLOATEST to keep system combined. Takes off modification 

° 
when o- is megabives 
Is O¢ — (7? Om positive? 


Store OL-A in O ww. 
(? 


Store Ge in Ow 


Is Oc = { negative? 

Oc - 7 Om-= 07 

Cc- 7? Ou = 07 

Is @e "ye ©Owmpositive? 

Is 9e “ Om positive? 

Modify FLOATEST to keep system combined. Takes off modification 
when ®,. is positive. 

Have 1 1/2 minutes of real time elapsed since the last print-out? 
Is Oc & 1.07 

set up to print when in the fourth quadrant. 

Cheek if in the foumth quadrant. 

Is the new ©O- greater than the previous C- 7 

Store thewnew time, C- and C, : 

EXIT (automatically set for desired jump-out). 

Print the last overshoot stored in LIMITBUF (symbol 888 ) 

Stop real time cilloek- 

Print position in phase trajectory where it was stopped by the 


real time clock (symbol 999). 
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88. 


89. 


90. 
91. 
92. 
99. 


99 - 
96. 
97. 
98. 


99 
100. 
101. 
aL Ze 
OD: 


104. 


Return to beginning to start new problem if new parameters 
remain to be evaluated. Unconditional stop if all parameters 
have been evaluated. 

Has the required number of overshoots been evaluated prior to 
starting the averaging routine? 

Add the @¢ just evaluated to the previous sum for averaging. 
Add one to the count prior to averaging. 

Clear time clock. 

Print time, 6, and Oe . 

Required number of total print-outs? 

Clear buffers. 

Obtain the average of G. . 

Peotracs 1.0. 

Print the average of O.~ 1.0 and how many overshoots were 
utilized to obtain the average. 

Is ©- less than 1.005 radians? 

Clear buffers. 

otart real time clock. 

EXIT (automatically set to desired jump-out). 

Has 40 seconds of real time elapsed since the last print-out? 
Print the position in the phase trajectory where it was stopped 
by the real time clock and ea. being less than 1.005 radians 


(symbol 1005). 
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Jump Switches: 

= - For LIMITCYC, allows the maximum overshoot print-out only. 
With J) down, a phase trajectory is printed for a maximun 
of le cycies. 

J - Set the number of print-outs depending on zeta. Change 
the number of overshoot values used to obtain the average 
value. 

J - Use the real time clock for an automatic recycle to a 
new parameter when have: 

a. 11/2 minutes maximum after the last print-out. 


= 


b. 40 seconds maximum and Oc is 1.005 radians or less. 
Stop Switches: 


é - Stops at the end of this parameter run or stops at the end 


of the manual print-out. 
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APPENDIX C 


PHASE TRAJECTORY PRINT OUT 


(Numbers indicated in tenths followed by power of 10) 
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